On Hr(curl,Ω) estimates for a Maxwell type system in convex domains

In bounded convex domains, the regularity of a vector field u with its $\mathrm{div}\mathbf{u}$, $\mathrm{curl}\mathbf{u}$ in $L^r$ space and the tangential component or the normal component of u over the boundary in $L^r$ space, is established for $1<r<\infty$. As an application, we derive an $H^r(\mathrm{curl},\mathrm{\Omega })$ estimate for solutions to a Maxwell type system with an inhomogeneous boundary condition in convex domains. In contrast to the well-posed region of r in the space $H^r(\mathrm{curl},\mathrm{\Omega })$ for the Maxwell type system in Lipschitz domains given by Kar and Sini (2016) [16], we extend the well-posed region to be optimal.     

Parametric error bounds for convex approximations of two-stage mixed-integer recourse models with a random second-stage cost vector

We consider two-stage recourse models with integer restrictions in the second stage. These models are typically non-convex and hence, hard to solve. There exist convex approximations of these models with accompanying error bounds. However, it is unclear how these error bounds depend on the distributions of the second-stage cost vector q. In this paper, we derive parametric error bounds whose dependence on the distribution of q is explicit: they scale linearly in the expected value of the ${?}_1$-norm of q.     

Hjelmslev quadrilateral central configurations

A Hjelmslev quadrilateral is a quadrilateral with two right angles at opposite vertices. Using mutual distances as coordinates, we show that any four-body central configuration forming a Hjelmslev quadrilateral must be a right kite configuration.     

Going Off the Grid: Iterative Model Selection for Biclustered Matrix Completion

We consider the problem of performing matrix completion with side information on row-by-row and column-by-column similarities. We build upon recent proposals for matrix estimation with smoothness constraints with respect to row and column graphs. We present a novel iterative procedure for directly minimizing an information criterion to select an appropriate amount of row and column smoothing, namely, to perform model selection. We also discuss how to exploit the special structure of the problem to scale up the estimation and model selection procedure via the Hutchinson estimator, combined with a stochastic Quasi-Newton approach. Supplementary material for this article is available online.     

Statistical and Geometrical Way of Model Selection for a Family of Subdivision Schemes

The objective of this article is to introduce a generalized algorithm to produce the m-point n-ary approximating subdivision schemes(for any integer m, n ≥ 2). The proposed algorithm has been derived from uniform B-spline blending functions. In particular, we study statistical and geometrical/traditional methods for the model selection and assessment for selecting a subdivision curve from the proposed family of schemes to model noisy and noisy free data. Moreover, we also discuss the deviation of subdivision curves generated by proposed family of schemes from convex polygonal curve. Furthermore, visual performances of the schemes have been presented to compare numerically the Gibbs oscillations with the existing family of schemes.     

Low rank matrix recovery from rank one measurements

We study the recovery of Hermitian low rank matrices $X\in {\mathbb{C}}^{n\times n}$ from undersampled measurements via nuclear norm minimization. We consider the particular scenario where the measurements are Frobenius inner products with random rank-one matrices of the form $a_j{a}_j^{?}$ for some measurement vectors $a_1,\dots ,a_m$, i.e., the measurements are given by $b_j=\mathrm{tr}(X{a}_j{a}_j^{?})$. The case where the matrix $X=x{x}^{?}$ to be recovered is of rank one reduces to the problem of phaseless estimation (from measurements $b_j=|〈x,a_j〉{|}^2$) via the PhaseLift approach, which has been introduced recently. We derive bounds for the number m of measurements that guarantee successful uniform recovery of Hermitian rank r matrices, either for the vectors $a_j$, $j=1,\dots ,m$, being chosen independently at random according to a standard Gaussian distribution, or $a_j$ being sampled independently from an (approximate) complex projective t-design with $t=4$. In the Gaussian case, we require $m\ge Crn$ measurements, while in the case of 4-designs we need $m\ge \mathit{Cr}n\log ?(n)$. Our results are uniform in the sense that one random choice of the measurement vectors $a_j$ guarantees recovery of all rank r-matrices simultaneously with high probability. Moreover, we prove robustness of recovery under perturbation of the measurements by noise. The result for approximate 4-designs generalizes and improves a recent bound on phase retrieval due to Gross, Krahmer and Kueng. In addition, it has applications in quantum state tomography. Our proofs employ the so-called bowling scheme which is based on recent ideas by Mendelson and Koltchinskii.     

Tensor Decomposition With Generalized Lasso Penalties

We present an approach for penalized tensor decomposition (PTD) that estimates smoothly varying latent factors in multiway data. This generalizes existing work on sparse tensor decomposition and penalized matrix decompositions, in a manner parallel to the generalized lasso for regression and smoothing problems. Our approach presents many nontrivial challenges at the intersection of modeling and computation, which are studied in detail. An efficient coordinate-wise optimization algorithm for PTD is presented, and its convergence properties are characterized. The method is applied both to simulated data and real data on flu hospitalizations in Texas and motion-capture data from video cameras. These results show that our penalized tensor decomposition can offer major improvements on existing methods for analyzing multiway data that exhibit smooth spatial or temporal features.     

Interval analysis method and convex models for impulsive response of structures with uncertain-but-bounded external loads

Two non-probabilistic, set-theoretical methods for determining the maximum and minimum impulsive responses of structures to uncertain-but-bounded impulses are presented. They are, respectively, based on the theories of interval mathematics and convex models. The uncertain-but-bounded impulses are assumed to be a convex set, hyper-rectangle or ellipsoid. For the two non-probabilistic methods, less prior information is required about the uncertain nature of impulses than the probabilistic model. Comparisons between the interval analysis method and the convex model, which are developed as an anti-optimization problem of finding the least favorable impulsive response and the most favorable impulsive response, are made through mathematical analyses and numerical calculations. The results of this study indicate that under the condition of the interval vector being determined from an ellipsoid containing the uncertain impulses, the width of the impulsive responses predicted by the interval analysis method is larger than that by the convex model; under the condition of the ellipsoid being determined from an interval vector containing the uncertain impulses, the width of the interval impulsive responses obtained by the interval analysis method is smaller than that by the convex model.The project supported by the National Outstanding Youth Science Foundation of China (10425208), the National Natural Science Foundation of China and Institute of Engineering Physics of China (10376002) The English text was polished by Keren Wang.     

Pathwise superhedging for time-dependent barrier options on càdlàg paths—Finite or infinite tradeable European,One-Touch,lookback or forward starting options

We establish pathwise duality using simple predictable trading strategies for the robust hedging problem associated with a barrier option whose payoff depends on the terminal level and the infimum of a càdlàg strictly positive stock price process, given tradeable European options at all strikes at a single maturity. The result allows for a significant dimension reduction in the computation of the superhedging cost, via an alternate lower-dimensional formulation of the primal problem as a convex optimization problem, which is qualitatively similar to the duality which was formally sketched using linear programming arguments in Duembgen and Rogers [10] for the case where we only consider continuous sample paths. The proof exploits a simplification of a classical result by Rogers (1993) which characterizes the attainable joint laws for the supremum and the drawdown of a uniformly integrable martingale (not necessarily continuous), combined with classical convex duality results from Rockefellar (1974) using paired spaces with compatible locally convex topologies and the Hahn–Banach theorem. We later adapt this result to include additional tradeable One-Touch options using the Kertz and Rösler (1990) condition. We also compute the superhedging cost when in the more realistic situation where there is only finite tradeable European options; for this case we obtain the full duality in the sense of quantile hedging as in Soner (2015), where the superhedge works with probability $1?\epsilon$ where $\epsilon$ can be arbitrarily small), and we obtain an upper bound for the true pathwise superhedging cost. In Section 5, we extend our analysis to include time-dependent barrier options using martingale coupling arguments, where we now have tradeable European options at both maturities at all strikes and tradeable forward starting options at all strikes. This set up is designed to approximate the more realistic situation where we have a finite number of tradeable Europeans at both maturities plus a finite number of tradeable forward starting options.¹     

An anisotropic area-preserving flow for convex plane curves

This paper will deal with an anisotropic area-preserving flow which keeps the convexity of the evolving curve and the limiting curve converges to a homothety of a symmetric smooth strictly convex plane curve.